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ABSTRACT 

Aims. While analyzing decades of very long baseline interferometry (VLBI) data, we detected the secular aberration drift of the extragalatic 
radio source proper motions caused by the rotation of the Solar System barycenter around the Galactic center. Our results agree with the 
predicted estimate to be 4-6 micro arcseconds per year (/ias/yr) towards a = 266° and 5 = -29°. In addition, we tried to detect the quadrupole 
systematics of the velocity field. 

Methods. The analysis method consisted of three steps. First, we analyzed geodetic and astrometric VLBI data to produce radio source 
coordinate time series. Second, we fitted proper motions of 555 sources with long observational histories over the period 1990-2010 to their 
respective coordinate time series. Finally, we fitted vector spherical harmonic components of degrees 1 and 2 to the proper motion field. 
Results. Within the error bars, the magnitude and the direction of the dipole component agree with predictions. The dipole vector has an 
amplitude of 6.4 ± 1.5 yuas/yr and is directed towards equatorial coordinates a = 263° and S = -20°. The quadrupole component has not been 
detected. The primordial gravitational wave density, integrated over a range of frequencies less than 10~ 9 Hz, has a limit of 0.0042 h~ 2 where h 
is the normalized Hubble constant is /f /(100 km/s). 
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1. Introduction 

Secular aberration drift is an apparent alteration in the ve- 
locity of distant objects caused by the acceleration of the 
Solar System barycenter in the Milky Way. The proper mo- 
tion field of extragalactic sources affected by the secular aber- 
ration drift should present a dipolar structure with amplitude 
4-6 microseconds of arc per year (yuas/yr) directed towards 
the Galactic center (see, e.g., Fanselow 1983, Bastian 1995, 
Eubanks et al. 1995, Sovers et al. 1998, Mignard 2002, 
Kovalevsky 2003, Kopeikin & Makarov 2006). 

Observations of extragalactic radio sources by very long 
baseline interferometry (VLBI) have been done since the late 
1970s by various agencies, mainly the United States Navy and 
the National Aeronautic and Space Administration (NASA), 
and have been coordinated since 1998 by the International 
VLBI Service for Geodesy and Astrometry (IVS, Schliiter & 
Behrend 2007). The main purpose of these observations is to 
monitor the Earth orientation and realize terrestrial and celes- 
tial reference frames. 

One of the cornerstones of the VLBI technique is the ac- 
curate realization of the celestial reference system using dis- 
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tant objects that supposedly have no detectable proper mo- 
tion (Feissel & Mignard 1998). The most recent realization 
of the International Celestial Reference Frame (ICRF2, Fey et 
al. 2009) provides absolute radio source coordinates for 3,414 
sources. For 1,448 ICRF2 sources observed in at least two ses- 
sions, the median formal error is 175 /ias. For the most ob- 
served radio sources, the inflated position errors are 40 yuas in 
each coordinate. 

This catalog did not include proper motions of sources be- 
cause the dominant proper motion for any source is the some- 
what temporary one associated with internal structure changes 
that produce apparent motions that are an order of magnitude 
greater than the aberrational drift (e.g., Fey et al. 1997). Such 
proper motions are not constant with time and are about an or- 
der of magnitude higher than the 4-6 yuas/yr proper motions 
expected from the secular aberration drift. However, we be- 
lieve that the secular aberration drift can be detected even with 
source position errors of ~ 100 yuas and apparent motion formal 
errors of 10-100 //as/yr (e.g., Feissel- Vernier 2003). 

VLBI measurement of the secular aberration drift is im- 
portant for both astrometry and astrophysics. Indeed, it would 
give an independent estimate of the acceleration vector ap- 
plied to the Solar System barycenter without refering to ob- 
jects within the Galaxy, so it can be used to constrain the mass 
of the central regions of the Milky Way. Moreover, the secu- 
lar aberration drift produces a slow deformation of the celestial 
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reference frame axes, and, if not corrected in geodetic VLBI 
sofware packages, could lead to contamination in estimates of 
other geodetic parameters, such as source coordinates and the 
Earth's nutation and precession (Titov 2010). 

Gwinn et al. (1997) have investigated the radio source ve- 
locity field obtained from the analysis of 16 years of VLBI data. 
Though the prime interest of their paper was to find an up- 
per limit for the mass-energy of the cosmological gravitational 
wave background from the quadrupolar structure of the veloc- 
ity field, their estimate of the Galactocentric acceleration is 
higher than expected by a factor of two with comparable stan- 
dard errors. MacMillan (2005) processed another seven years 
of data but did not detect this effect. Using the OCCAM geode- 
tic VLBI analysis software package, Titov (2009) reports a sta- 
tistically significant dipole harmonic, but its magnitude was far 
from the theoretical value. 

The approach used by MacMillan (2005) and Titov (2009) 
consisted of a direct estimation of the Galactocentric acceler- 
ation without intermediate estimation of the proper motions. 
Our approach differs from these two. We used a three-step 
procedure that includes (i) a production of source coordinate 
time series from analysis of VLBI delays, (ii) a least-square fit 
of proper motions to source coordinate time series, including 
editing large anomalous proper motions, and (iii) the fit of the 
Galactocentric acceleration to proper motions. This three-step 
analysis helped to remove outliers that significantly biased the 
results of the previous analyses. 

2. Aberration in proper motions 

Consider the Solar System barycenter moving on a quasi circu- 
lar orbit at distance R of the Galactic center with an accelera- 
tion 



d 2 R kM(R) 
a = —rr ; — "> 



dt 2 



(1) 



where k is the constant of gravitation, and M(R) the equivalent 
mass entering the problem. Using kM = u 2 ^ = V 2 R, one gets 
the acceleration a = dV/dt = -coVu, resulting in the aberration 
effect on the proper motion fi of distant bodies as seen from the 
Solar System barycenter (e.g., Kovalevsky 2003) 



u)V 
Afi = u, 

c 



(2) 



wherein c is the speed of light and u the unit vector pointing 
towards the Galactic center. 

Recent estimates of the Galactic parameters give a distance 
to the Galactic center of 8.4 + 0.6 kpc and a circular rotation 
speed of 254 + 16 km/s (Reid et al. 2009). The expected ac- 
celeration a = V 2 /R equals (2.5 + 0.5) x 10" 13 km/s 2 . This 
corresponds to an aberration of a distant body proper motion 
up to 5.0 + 1 .0 /ias/yr. 

Further accelerated motion of the whole Milky Way would 
produce an additional aberration. According to Doppler shift 
measurements of the Cosmic Background Explorer (COBE) 
satellite, which realizes a celestial reference frame based on 
the cosmic microwave background (CMB) radiation, the Local 
Group of galaxies moves at ~630 km/s (Kogut et al. 1993). 
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Fig. 1. The distribution of sources vs. the declination (Top), the 
proper motion fc 
nation (Bottom). 



proper motion formal errors cr = (o- 2 cosg + cr 2 ) 1 ^ 2 vs. the decli- 
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Fig. 2. The proper motion in right ascension vs. the right ascen- 
sion of 40 sources observed in more than 1,000 sessions. 
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Although considered as linear, this motion should correspond 
to a slowly accelerated motion on longer time scales, of orbit 
radius of several Mpc. It would result in an aberration that is 
far too small to be investigated here. 

For a distant body of equatorial coordinates (a, 6), © also 
reads as 

A// ff cos 6 = -d\ sin a + di cos a, (3) 
A//a = —d\ cos a sin 5 - di sin a sin 6 + di cos 6, (4) 

where the d; are the components of the acceleration vec- 
tor in unit of the proper motion, and which corresponds to 
degree 1 spheroidal (or electric) development of (see, e.g., 
Mathews 1981, Mignard & Morando 1990) 

i« = ZKA + < y fJ. (5) 
where d\ — of ,, dr — a? ,, d% — af n , and Yf and Yf are 

1 1,1' L 1,-1' J 1,0' l,m l,m 

the vector spherical harmonics of electric and magnetic types 
of degree I and order m. 

In addition to the aberration distortion, there may also be a 
small global rotation that can be described by the toroidal (or 
magnetic) harmonics of degree 1 : 

A/u a cos 5 - r\ cos a sin 5 + r2 sin a sin 5 - r; cos 6, (6) 
Aju,s = —r\ sin a + cos a, (7) 

where r, can be expressed in terms of vector spherical harmon- 
ics coefficients as r\ = af v = af_ v and = af Q . 

To investigate a possible quadrupolar anisotropy of the ve- 
locity field, we give the development of the degree 2 vector 
spherical harmonics (i.e., / = 2 in Eq. (O): 

A/u a cos 6 — — (a 2 2 e sin 2a — a^'^ cos 2a) cos 6 
+(«2 'j Re sin a - flj '[ Im cos a) sin 5 
+(af 2 Re sin 2a - af^ cos 2a) sin 8 cos 6 
+(a 2 w j Re sin a — a^f™ cos a) cos 25 
-#20 sin 6 cos 5, (8) 
A/us — -(flf? 6 cos ^ a + a 22 sm sm ^ cos 5 
-(«2 ' Re cos a + flf'j sin a) cos 25 

+ ^ fl 22 Re C0S ^ a + a 22 Im Sm C0S ^ 

-(a^f 6 cos a + a^ Im sin a) sin 5 

+«2 sin 5 cos 5. (9) 

3. Results and discussion 

3.7. Data processing 

We processed 5,030 sessions of the permanent geodetic and as- 
trometric VLBI program since 1979, totalling 7,285,312 group 
delay measurements at 8.4 GHz. Radio source coordinates 
were estimated once per session, together with Earth orien- 
tation parameters and station coordinates. The cut-off eleva- 
tion angle was set to 5°. A priori zenith delays were deter- 
mined from local pressure values (Saastamoinen 1972), which 



were then mapped to the elevation of the observation using 
the Vienna mapping functions (Bohm et al. 2006). Zenith 
wet delays were estimated as a continuous piecewise linear 
function at 30-min intervals. Troposphere gradients were es- 
timated as 8-hr east and north piecewise functions at all sta- 
tions except a set of 110 stations with short observational his- 
tories. Station heights were corrected for atmospheric pres- 
sure and oceanic tidal loading. The relevant loading quanti- 
ties were deduced from surface pressure grids from the U. S. 
NCEP/NCAR reanalysis project atmospheric global circula- 
tion model (Kalnay et al. 1996, Petrov & Boy 2004) and from 
the FES 2004 ocean tide model (Lyard et al. 2004). No-net 
rotation (NNR) and translation constraints per session were 
applied to the positions of all stations, excluding Fort Davis 
(Texas), Pie Town (New Mexico), Fairbanks (Alaska), and the 
TIGO antenna at Concepcion, Chile because of strong non lin- 
ear displacements (These two sites experienced post-seismic 
relaxation effects after large earthquakes on the Denali fault 
in 2003, and between Talca and Concepcion in early 2010). A 
priori precession and nutation comply with the IAU 2000/2006 
resolutions, which include the nutation model of Mathews et 
al. (2002), the improved precession model of Capitaine et 
al. (2003b), and the non rotating origin-based coordinate trans- 
formation between terrestrial and celestial coordinate systems 
(Capitaine et al. 2003a). Usually, an NNR constraint is used 
to fix the ICRS axes. However, Titov (2010) argues that ap- 
plication of a tight NNR constraint may wipe out all sys- 
temtic effects in the proper motion of reference radio sources. 
Therefore, we tied the celestial frame to the ICRF2 using a 
loose NNR constraint uniformly applied for each session. More 
details are discussed later in Section 3.2. The calculations used 
the Calc 10.0/Solve 2010.05.21 geodetic VLBI analysis soft- 
ware package, which was developed and maintained at NASA 
Goddard Space Flight Center, and were carried out at the Paris 
Observatory IVS Analysis Center (Gontier et al. 2008). 

Before 1990, the general deficiency of the VLBI networks, 
including the number of observed sources and observing an- 
tennas per session, makes the VLBI products less reliable (see, 
e.g., Gontier et al. 2001, Malkin 2004, Feissel-Vernier et al. 
2004, Lambert & Gontier 2009 who reports interesting statisti- 
cal results and remarks about the VLBI evolution over the past 
two decades). For this reason we removed data before 1990. A 
treatment of the full data base over 1979-2010 is nevertheless 
presented later for comparison. 

In the coordinate time series, data points resulting from 
fewer than three reliable observations within a session were 
removed, and outliers were eliminated so that the x 1 is rea- 
sonably close to unity. Then, proper motions were computed 
by weighted least-squares for time series containing at least ten 
points and longer than ten years. Weights were taken as the in- 
verse of the squared formal error. A set of 39 sources showing 
significant non linear positional variations due to large-scale 
variations in their structure (including 3C84, 3C273B, 3C279, 
3C345, 3C454.3, and 4C39.25) were isolated in the ICRF2 
work and treated in such a manner that they did not perturb 
the geodetic solutions (Fey et al. 2009). We removed these 39 
sources from our data set. 



4 



O. Titov et al.: VLBI measurement of the secular aberration drift 





DR 


DR1 


DR2 


DR3 


DR4 


No. sources 


555 


268 


257 


40 


515 


Dipole 












di 


-0.7 + 0.8 


-1.9 + 1.0 


-1.6 ± 0.9 


-3.4 ± 1.0 


4.2 ± 1.3 


d 2 


-5.9 ± 0.9 


-7.1 + 1.1 


-6.0 ± 1.0 


-6.2 ± 1.2 


-4.6 ± 1.3 




— 9 ? + i n 

— It 1 .u 


— 7. s + 1 9 
— J.J it 1 


— z..y it i.i 


— J.O It 1 .A 


14+17 

1 .H- It 1 . / 


Amplitude 


6.4 ± 1.5 


8.1 + 1.9 


6.9 ± 1.7 


8.0 ± 2.0 


6.4 + 2.5 


Rotation 












'•| 


-2.6 + 0.9 


—3.2 +1.2 


-2.7 + 1.1 


— 1 Q + 1 1 

i.y x i.j 


-4.3 + 1.4 


r 2 


0.4+1.0 


-0.6 + 1.2 


-0.3 ± 1.1 


2.5 ± 1.4 


-2.1 + 1.5 




0.8 + 0.7 


-0.3 ± 0.9 


0.0 ± 0.8 


-0.5 ± 0.9 


4.1 + 1.2 


Amplitude 


2.8 ± 1.5 


3.3 + 1.9 


2.7 ± 1.7 


3.2 ± 2.0 


6.3 + 2.4 


Direction of the acceleration vector 










Right Ascension (°) 


263 ± 1 1 


255 ± 11 


255 + 11 


241 ± 12 


312 ± 16 


Declination (°) 


-20 ± 12 


-25 ±11 


-25 + 11 


-28 + 12 


13 ±21 


Pre-fit wrms 


22.1 


19.7 


16.3 


9.2 


33.3 


Post-fit wrms 


21.7 


18.9 


15.7 


7.5 


32.7 


Reduced x 2 


1.9 


1.9 


1.8 


2.1 


1.9 



Table 1. The multipole coefficients (//as/yr) of the velocity field. Uncertainties are Icr. 



The final sample contains proper motions of 555 sources 
and is made available electronically. Figure [TJdisplays the dis- 
tribution of sources and proper motion formal errors in decli- 
nation. Near the polar areas, the number of sources decreases 
proportionally to the cosine of the declination. It also shows 
the nonuniformity of the sample and a lack of sources at de- 
clinations under -40°. Figure [2] displays /u a co&6 versus a for 
40 sources observed in more than 1,000 sessions (see Section 
3.2 for details). The apparent motions of these sources are esti- 
mated very accurately thanks to a large number of observations. 
A systematic in sin a of magnitude less than 10 /ias clearly 
shows up. Tiny underlying aberrational drift is indicative even 
for a limited number of well-observed radio sources. 

3.2. Dipole 

This section comprises our results of the dipole component es- 
timation. We start with a main solution including all 555 radio 
sources. Then we consider different subsets of radio sources to 
verify the robustness of the main solution. 

First, dipole and rotation coefficients were fitted by 
weighted least-squares following Eqs. (3)-(4) and (6)-(7) 
(Table [T] column DR). Reported errors are standard formal er- 
rors. The fit produces correlations of ~0.4 between di and r2 
and between ^2 and r\. Figure [3] displays both the proper mo- 
tions of the 555 sources and the estimated dipole component 
of the velocity field. Within error bars, the dipole amplitude 
and direction agree with predictions from measurements of the 
Galactic parameters. The corresponding centrifugal accelera- 
tion of the Solar System bary center is (3.2 + 0.7) x 10~ 13 km/s 2 . 



By fixing the distance of the Solar System to the Galactic cen- 
ter to 8.4 kpc, the DR solution yields a circular rotation speed 
of 276 km/s. 

To test the robustness of our estimates, we fitted the param- 
eters to various subsamples of sources, including sources with 
structure indices (Chariot 1990, Fey & Chariot 2000) lower 
than 2.5 (DR1), thus keeping only the most compact sources 
only. Since structure indices were not available for all sources, 
the number of sources in the sample is considerably reduced. 
Nevertheless, the amplitude and the orientation of the dipole 
do not change significantly with respect to DR. Another test 
consisted of only keeping ICRF2 defining sources (DR2). Two 
hundred fifty seven of them have a sufficient observational his- 
tory to pass our selection scheme successfully. Again, it does 
not change the amplitude and the orientation of the dipole sig- 
nificantly. 

In all the above solutions, the rotation is about 5 /ias/yr and 
is consistent with the ICRF2 error. If the rotation is fixed to 
zero, the estimated dipole is of 7.1 + 1.5 /zas/yr toward (a - 
263 + 9°, 5 = -17+ 11°). 

Owing to the uneven number of observations for the 555 
radio sources, their proper motion standard errors vary by three 
orders of magnitude: from 5 to 5000 /ias/yr. It is necessary 
to check whether the estimates of the dipole component are 
dominated by a small number of radio sources with small er- 
rors on proper motion. Another fit was therefore done to the 
40 sources observed in more than 1,000 sessions and used to 
produce Fig. [2] (DR3). These sources have formal errors of 
proper motions less than 15 //as. DR3 has the best weighted 
rms among all the solutions in Table [TJ 
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NNR cr Amplitude a 6 



Yes 


2 as 


6.4 


± 1.5 


263 ± 11 


-20 ± 


12 


Yes 


2 mas 


5.2 


±0.7 


237 ±7 


22 


+ 8 


No 


2 as 


6.5 


± 1.5 


263 ± 11 


-23 ± 


12 


No 


2 mas 


3.3 


±0.5 


242 ±7 


22 


+ 7 



Table 2. The dipole amplitude (ji/as/yr) and orientation (°) ob- 
tained with different constraints with uncertainties of lcr. 



The last fit (DR4) was done to the remaining 515 sources 
observed in less than 1,000 sessions. DR3 and DR4 deviate 
from the main solution DR in opposite directions. Though 
the orientation of the dipole is mainly constrained by the 40 
sources with the longest observational history, the magnitude 
of the dipole is closer to the theoretical values when the 515 
sources are used. 

We also made the adjustment without removing data before 
1990. It yielded a dipole amplitude of 6.4 ±1.5 //as/yr, directed 
towards (a = 262 + 10°, 5 = -15 ± 12°). Though the amplitude 
is comparable to the one obtained over 1990-2010, the direc- 
tion appeared changed by almost 5° in declination. It is likely 
due to a lack of sources observed in the southern hemisphere 
during the first decade of VLBI. 

Overall, the estimate of the dipole effect obtained from the 
DR solution (6.4 ± 1.5 //as/yr; a = 263 ± 1 1°, 5 = 20 ± 12°) 
is considered confident and could be recommended for futher 
applications. 

3.3. NNR constraint 

Before going further, we address the problem of the constraint 
that was applied to the observation during the first step of the 
VLBI data reduction (i.e., production of source coordinate time 
series). Table|2]displays the amplitude and the orientation of the 
dipole fitted to time series obtained with different constraints. 
The first line is the DR solution for which the cr of the NNR 
constraint is 2 as. Decreasing the cr to a value close to the milli 
arcsecond or lower makes the dipole move off the Galactic cen- 
ter. 

The last two lines of Table [2] mention solutions without 
NNR. Applying no constraint at all would cause degeneracy 
of the normal equation matrix for daily parameters. Therefore, 
in this solution, each source is tied to the ICRF2 position using 
a loose constraint of 2 as. We can see that such a constraint is 
equivalent to the NNR applied with the same cr. 

3.4. Quadrupole 

In a second fit, we adjusted all parameters of Eqs. (3)-(4) and 
(6)-(9) (Table[3]l. Maximum correlations of 0.5 showed up be- 
tween d\ and a^j Re . Including of the quadrupole harmonics did 
not change the dipole and rotational effects. The amplitude of 
the aberrational effect increased by 0.7 yuas/yr, and the right as- 
cension of the vector moved off the theoretical value, but the 
declination of the vector direction almost coincided now with 
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Fig. 3. {Top) The proper motions of the 555 sources, and 
(Bottom) the estimated dipole component of the velocity field. 
The green line represents the equator of the Milky Way, whose 
center is indicated by the blue marker. 

the theoretical prediction. Fixing the dipole and rotation param- 
eters to the DR values and adjusting the quadrupole parameters 
only leads to 6.4 + 3.4 //as. 

The quadrupole component may come either from the 
Hubble constant anisotropy or the primordial gravitational 
waves (Kristian & Sachs 1966, Pyne et al. 1996, Gwinn et 
al. 1997, Jaffe 2004). Though the Hubble constant anisotropy 
only affects af and a^ e (Titov 2009), the primordial grav- 
itational waves affect all the quadrupole terms in Table [3] 
Whereas the resultant amplitude is less than 3cr standard error, 
the only marginally statistically significant component a 2 \ e = 
4.1 + 1.3 //as/yr can be converted into the energy density of grav- 
itational waves Qgw- Using Eq. (1 1) of Gwinn et al. (1997) for 
this component, one gets Qgw = 0.0042 + 0.0004 hr 1 , where 
h = //o/( 100 km/s) is the normalized Hubble constant. The 
squared proper motion of quasars at cosmologic distances are 
sensitive to gravitational waves of long wavelength, even com- 
parable to the scale of the Universe, and proportional to Qgw 
on a wide range of frequencies, from the inverse of the pe- 
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No. sources 


555 


Dipole 




d, 


0.7 ± 0.9 


di 


-6.2 ± 1.0 


u 3 


-3.3 ± 1.0 


Amplitude 


7.1 ± 1.7 


Rotation 




n 


-2.4 ± 1.0 




0.4 ± 1.1 




0.8 ± 0.7 


Amplitude 


2.6 ± 1.7 



30 r 



• Dipole 

A Quadrupole 
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= 1.2 








a M 
"2,2 


-0.6 = 


= 1.3 


2.4 


+ 


1.3 


"2.1 


2.0 = 


= 1.1 


-0.6 


+ 


1.1 


"2,0 


0.7 = 


= 0.9 








Amplitude 


6.4 = 


= 3.6 









Direction of the acceleration vector 


Right Ascension (°) 


277 ± 12 


Declination (°) 


-28 ± 10 


Pre-fit wrms 


22.1 


Post-fit wrms 


21.6 


Reduced x 1 


1.9 



Table 3. The quadrupole coefficients (/^as/yr) of the velocity 
field with uncertainties od lcr. 



riod of observations to Hubble time. Therefore, the value es- 
timated above may indicate the upper limit of the gravitational 
waves density integrated over a range of frequencies less than 

irr 9 Hz. 

3.5. Dependence on the redshift 

We also looked at the dipole and quadrupole amplitudes as 
functions of the redshift, available for 488 sources (Titov & 
Malkin 2009). Each estimate incorporated a subsample con- 
taining 122 sources whose redshifts are between the values 
indicated by the x-axis ticks of Fig. [4] The mean declination 
of each subsample is between 7.2° and 11.7°. The dipole and 
quadrupole amplitudes do not present significant dependence 
on z. 

4. Conclusion 

This study showed that VLBI has now accumulated accurate 
enough data to detect the Galactocentric acceleration through 
its effect on distant radio source positions. It turns out that the 
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Fig. 4. The dipole (red disks) and quadrupole (blue triangles) 
amplitudes as functions of the redshift. Uncertainties are lcr. 

current definition of the celestial reference frame as epochless 
and based on the assumption that quasars have no detectable 
proper motions should be mitigated. In the future, VLBI real- 
izations of the celestial reference system should correct source 
coordinates for this effect, possibly by providing source posi- 
tions, together with a corrective formula. 

The European optical astrometry mission Gaia (Perryman 
et al. 2001), scheduled for 2012, should be able to determine 
the components of the acceleration vector with a relative pre- 
cision of 10% (Mignard 2002). To improve the VLBI deter- 
mination of the Galactocentric acceleration and to confirm the 
significance of the quadrupole systematics, more proper mo- 
tions of extragalactic radio sources need to be measured over 
the next decade. Concentrating on sources showing a high po- 
sitional stability and having a low structure index would re- 
duce unwanted effects of intrinsic motion caused by the rela- 
tivistic jets and other modification of the source structure. In 
addition, it is necessary to run a dedicated program to measure 
the redshift of the reference radio sources using large optical 
facilities, especially in the southern hemisphere (Maslennikov 
et al. 2010). 
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